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We study the regularization dependence on meson properties and the phase diagram of quark 
matter by using the two flavor Nambu-Jona-Lasinio model. We find that the meson properties and 
the phase structure do not show drastically difference depending the regularization procedures. We 
also find that the location or the existence of the critical end point highly depends on the regular¬ 
ization methods and the model parameters. Then we think that regularization and parameters are 
carefully considered when one investigates the QCD critical end point in the effective model studies. 
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I. INTRODUCTION 

The phase structure of quark matter on finite temperature and density has actively been studied for decades [jj. 
Under usual condition, meaning low temperature and density, quarks are confined inside hadrons and they never be 
able to observed as a single particle. On the other hand, due to the nature of the asymptotic freedom |2(, quarks 
and gluons can be free from the confinement at high temperature and density, because the coupling strength becomes 
weak at high energy. It is, therefore, expected that quark matter undergoes the confined/deconfined phase transition 
at some temperature and density. This is important subject both in theoretical and experimental studies since it 
crucially relates to the quark matter properties at relativistically high energy collisions and extremely dense stellar 
objects such as neutron stars. 

The first principle for quarks and gluons is Quantum Chromodynamics (QCD) which is a non-Abelian gauge field 
theory for fermions. Our goals is to evaluate the phase structure based on this first principle QCD, however, it is 
difficult to extract theoretical predictions due to the nature of complicated strongly interacting system. One of the 
most reliable approaches is to use the discretised version of QCD called the Lattice QCD (LQCD) in which theoretical 
calculation is performed on the discrete spacetime Although the LQCD works well at finite temperature T for 
small chemical potential /i ~ 0, there is the technical difficulty called the “sign problem” when one tries to investigate 
the system at intermediate chemical potential. There effective models maybe nicely adopted because some models 
can consistently treat the system at finite temperature and chemical potential. 

For the sake of evaluating the phase structure of quark matter at finite temperature and chemical potential, we will 
employ the Nambu Jona-Lasinio (NJL) model 4 which is the most frequently used one in this context (there are a 
lot of nice review papers on the model, see, e.g. (Bj-O].) The model is constructed by incorporating the four point quark 
interaction into the model Lagrangian, so it is not renormalizable due to this higher dimensional operator. Therefore, 
the physical predictions of the model inevitably depend on the regularization procedure and the model parameters 
chosen. The resulting phase diagram on the T — n plane is as well affected by the parameters and regularization 
prescriptions. So it is an important issue to study whether the phase structure obtained in one regularization method 
is consistent with the ones from different regularization methods. 

In this paper, we are going to study the phase structure of quark matter in the NJL model with various regularization 
ways, which are three dimensional (3D) momentum cutoff, four dimensional (4D) momentum cutoff, Pauli-Villars (PV) 
regularization, proper-time (PT) regularization, and the dimensional regularization (DR). The 3D cutoff scheme is the 
most popular method in this model and a lot of works have been done in this way. The 4D cutoff method preserves 
the Lorentz symmetry in which space and time are treated on equal footing. The Pauli-Villars regularization is based 
on the subtraction of the amplitude conside ring the virtually heavy particle to suppress the unphysical high energy 
contribution coming from loop integrals fl(J4l2|. The proper-time regularization makes integrals finite through the 
exponentially dumping factor [111 . 1131 ]. The dimensional regularization analytically continues the spacetime dimension 
in the loop integrals to a non-integer value, then try to obtain finite contribution from the inte gra ls M- Beside from 
the frequently used 3D cutoff way, there have been a lot of works by using the 4D [(J, @! , PV f6l. PT j(| ITsl [2(i} , 

and DR [2714321 ] . The physical consequences depend on the regularization fUJ. 

This paper is organized as follows; Section [TT] introduces the model Lagrangian, and show the model treatment on 
the meson properties and the explicit formalism at finite temperature and chemical potential. In Sec. |ml we present 
various regularization procedures, 3D, 4D, PV, PT and DR prescriptions with explicit equations. We then perform 
the parameter fitting in Sec. IIVI In Sec. |Vj the numerical results of the meson properties are shown. We then draw 
the phase diagrams with several parameter sets using various regularization methods in Sec. IVII We also study the 
phase diagram with the parameters fixed under the condition with the same constituent quark mass in Sec. IVIII 
In Sec. IVIII1 we give the discussions on the obtained results. Finally, we write the concluding remarks in Sec. IIXI 
Several detailed calculations are shown in Appendix. 


II. TWO FLAVOR NJL MODEL 

In this paper we consider two light quarks with equal mass. The model has SUl( 2) <E> SUr( 2) flavor symmetry at 
the massless limit, m —> 0. 


A. The Lagrangian and gap equation 

The Lagrangian of the two flavor NJL model is given by 

C = 4>(ifi - rh)ip + G [(tpip) 2 + (^*75T a ^) 2 ] , 


( 1 ) 
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where to is the diagonal mass matrix to = diag (m u ,md) and G is the effective coupling strength of the four point 
interaction. We set = m u in this paper. The application of the mean-field approximation 

( 2 ) 

leads the following mean-field Lagrangian 

2 

£ = $(i ft - m*)ip - (3) 

with the constituent mass to* = m u + cr. The flavor symmetry is broken down, S'L r i(2) ® SUr(2) —► SUl+r(2), by 
non-vanishing current quark mass, m u , and dynamically generated a. Thanks to the simple form of the Lagrangian, 
one can easily evaluate the effective potential, V G ff = — In Z/V where Z is the partition function 


Z= T>[ijj] exp 


i / d 4 xC 


and V is the volume of the system. After some algebra, we see 

~ 2 r cl A k 


= 4G “ 


*( 2tt) 4 


lndet(# — to*). 


The detailed derivation of the effective potential is presented in [9(. 

The gap equation is obtained through the extreme condition of the potential with respect to a , namely, 


dVeff 

da 


= 0. 


This condition leads the following gap equation 


a = 21V fG ■ 


with the number of flavors Nf and 


itrS(m*) = — tr J — 


d 4 k 


(27r) 4 ji — to* + is 1 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


where trace takes the spinor and color indices. This is the key equation in the model because it determines the values 
of the chiral condensate (ipip) and the constituent quark mass to* . 


B. Meson properties 


The properties of the pion and sigma meson can be studied based on the model with the determined chiral conden¬ 
sate. The interacting Lagrangian of the pion and quarks is written by 

C-tc qq = i9nqqi^/5T ' (9) 

where r,; are 2x2 matrices in the flavor space and ir z represent the pion fields. The explicit expression is r • 7r = 
t_7t - + t_|_7t + + To7T°, with t± = (tj ± T2) / y/2 and ro = t 3 where Tj are the Pauli matrices. 

By applying the random phase approximation, we can write the pion propagator as the summation of the geometrical 
series of the one-loop diagram, which gives 


A AP 2 ) = 

where IL 7r is the following quark loop contribution 


yirqq 


2 G 


p 2 — m 2 1 — 2GIT r (p 2 ) 


n ”(p z ) = -2 


f d 4 k 

i( 2n) 4 


tr [j 5 S(k)"/ 5 S{k -p)] 


( 10 ) 


( 11 ) 


with the quark propagator 


S(k) 


1 

ji — to* + ie 


( 12 ) 
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The explicit derivation of Eq. m is discussed in the review paper [g]. The pion mass is calculated at the pole position 
of the propagator, so the condition reads 


[1 - 2 GW(p 2 )\ 



(13) 


It should be noted that the residue at the pole p 2 
have the relation 
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2 

irqq 


m 2 coincides with the square of the coupling strength g 2 qq so we 



(14) 


In the similar manner, the sigma meson mass is evaluated at the pole of the propagator, 

A (p 2 ) = 9 ° qq _ — _ 

ayP ’ i?-m* _ l-2GII CT (p 2 )’ 

with 

/ d 4 /r 

—tr[S(k)S(k-p)]. 

Therefore the condition which determines the sigma meson mass becomes 

[1-2GITV)] U^=0. 

The pion decay constant is calculated from the following equation 

iSijP^fn = ( 0 | 


(15) 


(16) 

(17) 

(18) 


The explicit form becomes 

1 r d 4 k 

P p U =2 J tr b P ^9^qqS{k)^S{k - p )\. (19) 

Thus the equations Eqs. m, on and (11911 are the ones which determine the pion mass, sigma mass, and the pion 
decay constant. 


C. Explicit formalism at finite temperature 


Since our purpose is to study the phase structure on temperature T and chemical potential /i, we need to extend 
the equations to finite temperature. According to the imaginary time formalism, the integral region of the temporal 
component becomes finite due to the periodic or anti-periodic condition of fields as 


Z = 




( 20 ) 


where r is imaginary time and /3 is the inverse temperature 1 /T. Consequently, continuous integral in the temporal 
direction is replaced by the following discrete summation, 


d 4 fc 

i{2nY 


■F(k 0 ,k)^T J2 


n =—oo 


d 3 k 


F(iu n + p, k), 


( 21 ) 


where uj n = 2nirT or (2n + 1)7 tT depending on the statistical property of field, i.e., for bosons or fermions, and the 
chemical potential seen in Eq. (1201) appears in the way iu n + p. In this paper, we only treat fermionic quark loop 
contributions then u n = (2n + 1)7 tT is always the case. 

With the help of the formalism Eg. (1711) . we see that the gap equation at finite temperature becomes 


tr = 2N f G ■ [trS° + trS T ], 


tr S° = -N c m* , 


’ d 3 fc 1 
(2 tt) 3 2E’ 

d 3 k 1 


E/(^) 


( 22 ) 

(23) 


tr S T = N c rri 


(2tt) 3 2 E 


(24) 
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where N c is the number of colors, E = \/k 2 + to* 2 , E ± = E± p and f(E) = 1/(1 + e^ E ). It is important to note that 
the contributions can be expressed by the summation of the T independent part (trS 10 ) and T dependent part (tr,S' T ). 
This characteristic is general if one takes the infinite number of the frequency summation in finite temperature field 
theory and crucial when we apply the regularization procedures to the appearing integrals. 

Since the gap equation is derived by differentiating the effective potential with respect to a 1 then the effective 
potential can be obtained by integrating the gap equation (see, for example, HI), 

V(tr) = ^ - iV/ f da' itrS(m u +a'). (25) 

4(j Jo 

where we have dropped the suffix in V c g and just written V for notational simplicity. Thereafter the effective potential 
at finite temperature V = V a + V° + V T is evaluated as 


V a = 


a 


2 


4 G’ 


V° = 



d 3 fc 

(2/bj 3 


E , 


V T = -2 N f N c T 



(26) 

(27) 

(28) 


It is important to note that, if we apply some regularizations, the results between the direct calculation from Eq. © 
and the one after integrating the gap equation may become different, because regularization essentially means the 
subtraction and there are several ways of subtractions. Therefore, in this paper, we persistently use the latter way 
shown in Eq. (E5l) so that the model treatment becomes consistent. It should be noticed that the finite temperature 
correction, V T , contains no divergent integral. A finite result can be obtained for the finite temperature correction 
without applying any regularizations. 

Next, we carry on the integral in the meson properties; the one loop contribution can be written as 


nV) = -^ + P 2 /(p 2 ), 

m 


with 


I(p 2 ) = W c J - 


d 4 fc 


(27t) 4 (k 2 — m* 2 )[{k — p) 2 — m* 2 ] 


(29) 

(30) 


Since trS 1 is already evaluated above, the remaining task is to calculate /(= 1° + I T ), and it becomes 

(31) 

(32) 

Similarly, the one-loop diagram of the scalar channel can be written as 

rrV) = + (p 2 - 4m* 2 ) I{p 2 ). (33) 


\ P 2 ) = m c J 
r ( P 2 ) = -m c J 


i 


(27r) 3 E(4E 2 - p 2 ) ’ 
d 3 k £± f(E±) 
(2tt) 3 E{4E 2 -p 2 )' 


We now have already evaluated all the ingredients of H a above in Eqs. m, <ED and so we do not need 

further calculations. 

Finally, let us derive the equation for the pion decay constant. After a bit of algebra we obtain the relation, 


Here we evaluate at p 2 = 0 following [g]. 


/?r — 7(0). 


(34) 


III. REGULARIZATION PROCEDURES 

Since the integrals obtained in the previous section include infinities, we need to apply some regularization so that 
the model leads finite quantities. As mentioned in the introduction, the model is not renormalizable, then the model 
predictions inevitably depend on regularization procedures chosen. Here, we shall present possible regularization 
methods in this section. 
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A. Three dimensional cutoff scheme 


The idea of the three dimensional (3D) cutoff is simple; one drops high frequency mode by introducing the cutoff 
scale A 3D into the integrals. We work in the 3-dimensional polar coordinate system and cut the radial coordinate as 


d 4 fc r dk 0 f AsD k 2 dk 

(2tt) 4 J 2tt J 0 (2?r) 3 


dD 3 . 


By performing the integrals, we have for the gap equation a = 2NfGtrS, 

N c m* 


f r S ' 0 - 

“ °3D — 


27T 2 


- m*“ In A3D + VA^ + m^ 


— — 


N c m* f A3D 


d k 


E 


E/( E± ) 


The effective potential can also be calculated as 


v o M NcNf 
V3d((T) - 


A 3 d \j- A 3D + m* 2 (2A 3D + m* 2 ) - ro* 4 In A3D ^ A 3D + m_ 


^3d( (J ) — ~ 


N c NfT [ Asd 


d k k l 




The quark loop integral in the meson properties I(p 2 ) reads 

J 0 D = ^/ A 3 D d fc * 


7T 


fO 


js,-^ r-d* 


E(4E 2 - p 2 ) ’ 

k 2 


7T 


to 


E{4E 2 - p 2 ) 


E/(^) 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


Note that the integral diverges around 4 E 2 ~ p 2 , and we apply the principal integral to avoid this divergence HU- It 
may be worth mentioning that the integral can be performed analytically when p 2 = 0 for T = 0, then one has for 
the pion decay constant, 


/7r3D — 


N c m 

2n 2 


*2 


A 


3D 


- In ■ 


A; 


3D 




3D 


m 


*2 


V A 3D + m * 2 

We thus obtain the required quantities in evaluating the phase diagram and meson properties. 


(42) 


B. Four dimensional cutoff scheme 


In the four dimensional (4D) cutoff regularization scheme, we introduce the cutoff scale A 44 ) in the Euclidean space 
after performing the Wick rotation, 


/ 


d 4 k E 


f AiD k%dk E 

Jo (27T) 4 


d£l±. 


(43) 


This is well known four dimensional cutoff method for T — 0 case. As the natural extension to finite temperature, we 
introduce the cutoff scale by 


f d 4 k E y- /VAJn-o* k 2 dk 

J J^Jo W 


dfl 3 , 


(44) 


where L 4 is the maximum integer which does not exceed A^ / (2nT) — 1/2. 

In the 4D cutoff way, it is difficult to divide the contribution into the temperature independent and dependent 
parts, since there is also cutoff in the frequency summation. 
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The explicit form of tr S and the effective potential become 


trSho 


V4D (cr) 


N c m*T 

2tt 2 

N c NfT 

47T 2 


E 

n=— L4 — 1 J u 

L 4 

E 

n=— La — 1 * 


r*\/^4D ^ 


dkk z 


(LUn) 2 +E 2 ’ 




dkk ln((w„) +E ), 


(45) 

(46) 


where uj n = oj n — ip. 

For T = 0, the integral can be performed analytically by using the Feynman parameter method, 


tr <? u - 

ud 4 d — 


N c m* 


^ > 4D( cr ) — ~ 


7T Z 

N c N f 
8 t r 2 


A^p) TO 


*2 


in 


^4D 


m 


*2 


A 


4D 


*2 

m — m 


*4 


m 

A 2 


*2 


In- 


4D 


^*2 


A 4 d ln(A 2 


4D 


m 


* 2 \ 


(47) 

(48) 


One should give the special attention in calculating / 4 d(p 2 ), because the integral includes divergence to be cured as 
seen in the 3D cutoff case. The analytic expression of / 4 D (p 2 ) will be given in appendix [A] 

Again we show the explicit form for the pion decay constant at T = 0, 


fir 4D — 


N c m 
47r 2 


*2 


A 2 

]1 4D 


^4D 


1 *2 


In / A 4° 


t *2 \ 1 


*2 


(49) 


C. Pauli-Villars regularization 


In this regularization, the divergences from loop integrals are subtracted by introducing virtually heavy particles as 


1 1 



E 


at 

k 2 - A 2 ' 


(50) 


This manipulation induces virtual frictional force so that the contribution from unphysical high frequency mode is 
suppressed. 

In evaluation the gap equation, we apply the following subtraction 


1 Ql 02 

p 2 — to * 2 p 2 — A 2 p 2 — A 2 ’ 

where 

to * 2 — A| A 2 — TO* 2 

ai = A 2 — A 2 ° 2 = A 2 — A 2 


By setting the cutoff scales A 4 = A 2 = Apy after the subtraction, we have 


tl’S'py — 


N c m* 

47T 2 


Apy — to * 2 + to * 2 In 


,*2 


A * 2 

yl PV 


tiSpy = —2 N c m* 


d 3 p 

( 2 ^ 


f( E m) _ L , Apy — TO. 


E„ 


*2 


2p 2 


f(E 


±M 


E a 


(51) 


(52) 


(53) 

(54) 


where E m = y/k 2 + to * 2 and E\ = \Jk 2 + Ap V . By integrating the above equation, we obtain the effective potential 

N c N f 


V° — —- 

1 /py — 


V PV — — 


8 tt 2 
N c N f T 


A 2 *2 3 *4 . 1 *4 1 TO 

ApyTO - -TO + -TO In ^2 


(55) 


/ d/c /c 2 ln(l 




PVJ 

TO 


)-(4fc 2 + 2A 2 - m* 2 ) 




Ea 


(56) 
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Since the divergence coming from the integral I{p 2 ) is order of log, one subtraction is enough to make it finite, so we 
get 


T° — 

1 PV — 


2 N r 


7T 


[I - _ 
1 PV — 


2 N c 

7T 2 


dkk 2 


dkk‘ 


E m {dE 2 m -p 2 ) E A (dE\-p 2 ) 

E± f{E£) E± f{E±) 

E m {dE 2 m -p 2 ) E A {dE\-p 2 ) 


The pion decay constant at T = 0 becomes 

c2 N c m * 2 
- 4?r 2 


-1 + 


A 2 

2V pv 


Apy to 


*2 


A 2 

In^Y 


*2 


(57) 

(58) 


(59) 


D. Proper-time regularization 

The basic idea of the proper-time regularization is based on the following manipulation of the Gamma function, 

dr r n ~ 1 e~ Ar , (60) 


1 1 

A n T[n] 


' * /^p-i 


where the lower cut l/Ap T induces the dumping factor into the original propagator since, for example with n = 1, 


X E 


dr e At = 


1/A| n 


k% 


o C k E +m )/A pT 


(61) 


Therefore in this regularization high frequency contribution is dumped by the factor e fc £/ A PT, so the original divergent 
integral turns out to be finite. For A contains a imaginary part, Eq. is modified as, 


I jn /-co 

ffn] J o dT r " _le_lAT ’ ( Im (A) < 0,Re(n) > 0). 


Under this procedure, the integral of tr S in the gap equation becomes 


o _ N c m* 
PT “ ~AE r 


[Ap X e -m * 2/A P T +m* 2 Ei(-m* 2 / Ap T ) 


tr5pT — 


N c m*T 


E 


dr 


2 - 3/2 t'oho r 3/2 


; -i7r/4 e -i{(o; n ) 2 +m‘ 2 }r + cc ( 


(62) 

(63) 

(64) 


where Ei(—x ) is the exponential-integral function. For m * 2 -C Ap T , Eq. (16311 is expanded as 


tr9 o „ N c m* 

tiOp X _ 


Ap T — to* + to* I In —j-h 7 e — 

V Apt 

We rotate the contour of the integration in Eq. (1M1) to the imaginary axis of r 
trace becomes 


TO 


*2 


2A| t 


N c m*T “ 

trOpT = -777;— / 


-3/2 


n —0 - L / il PT 


dr 

73/2 


cos(2w n ^tr)e 


— (iai 2 — /x 2 +m* 2 )r 


(65) 

. For Wg — p , 2 + to* 2 > 0, the 

, ( 66 ) 


and for ujq — p 2 + to* 2 < 0, 


trSpx = 


N c m*T 

tt3/2 


E 

n>[iV] " ‘('‘PT 


dr 

' — 3 J 2 cos( 2 u) n pT)e 

i/A.1 t*/ 2 


( 2 2 1 *2\ 

— \^n —M ~\~TYh ) 


IN] 

n =0 ^* /i /- /v -pT 


00 dr . . 

sm(2 u} n pT)e 

i-tv/2 


fi 2 +m* 2 )T 


- Apt Re ( e™ /A f d6 e~ i6/2 exp [-/{(w”) 2 + m* 2 }e lf> /Al T \ 

V d-^/2 


( 67 ) 
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where N = {yj[i 2 — to* 2 /( 7 rT) — l}/2 . Similarly, the effective potential can be calculated through 


Vpt(c)= 

For ujq — fj, 2 + to* 2 > 0, one has 


N c NfT y , j dr r -3i7T/4 ) 2 +m* 2 |T . 


n=0 ' 71 / A pt 


( 68 ) 


(69) 


and for Wg — /j, 2 + to* 2 < 0, 


,, , x JV C IV>T 

V PT(a) - 


E 


dr 


n>[JV] J 1 / A pt T 
dr 

n—0 I-^/Apt 


.5/2 


cos(2 uj n fj,T)e M +m ^ 


1I.C 


+ A 3 Re fe" l7r/4 J ' dO e~ 3l6/2 exp [-i{(w“) 2 + m* 2 }e ifl /A£ T ] 


IpT can also be calculated by 


At t TfL r 1 / 2 

/pt(p 2 ) = -^lE / da f 

n n=0 J ° J +° 


dr 

2 


e -3i7r/4 e -i{(o, n ) 2 +A}r +cc _ 


(70) 


(71) 


where a is the Feynman integration parameter and A = to* 2 — p 2 /4 + a 2 p 2 . Then the integral can be written 


9 At T 7 1 / 2 , oc . 

/pT (P 2 ) = T372- [ da E 


n—0 


/*°° 2 2 
9(W n (a)) - T75 -cos(2w n /xr)e _(u; " _M +A)t 

9 1/ApT ^ 1/2 


dr 


+9(-W n (a)) < / —-7T- sin(2w n ^T)e 

1-A/a 2 t t 1 / 2 


(w 2 -/x 2 +A)t 


- i Re U“ i7r / 4 J 7 d<9 e ie / 2 exp [-*{(w") 2 + A}e ie /Aj3 T ] 


(72) 


where W„(q) = w 2 + to* 2 + (a 2 — l/4)p 2 — /r 2 . 

The pion decay constant at T = 0 reads the following simple form, 


N, —* 2 


fZ P T = - I -^Ei(-m* 2 /A 2 PT ). 


For to* 2 -C Ap T , we have 


fl pt — 


N r m 


,*2 


47r 2 




„*2 


A 2 

i\_PT 


2 /A 2 t ). 

(73) 

+ ln A (A. 

m z J 

(74) 


E. Dimensional regularization 


In the dimensional regularization method, we obtain finite quantities through analytically continuing the dimension 
in the loop integral to a non-integer value, D, as 


d 4 fc 
(2^^ M ° 


4—D 


d D k 
(27r) 


D ’ 


(75) 


where the scale parameter Mq is inserted so as to adjust the mass dimension of physical quantities. The method is 
well known since it preserves most of symmetries. Note that this result is the same as the result obtained from the 
proper-time integral (0 < r < 00 ) and expressed by the poles of the Gamma function. 
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The trace tr S in the gap equation reads 


t cO _-N c Mt°m * 
wd dr — 


(:2 n) D /2 


trS'op, = — A D m* / dfc- 


r (l - y) 


'„D -2 


where 


The effective potential becomes 


N C 2 2 ~ D I 2 M^~ D 
D ~ 7r(£>-i)/2r((D - l)/2)' 


V 0 _ N c N fM% d / *2\D/2 

Vdr " 2(2 tt)W 2 H 2; (m } ’ 

V£ r = —A D N f T f dfc fc 15 " 2 ^ In [l + e _/3£; 

In the similar manner, the integral /(p 2 ) appearing in the meson propagator is calculated as 

' k D ~ 2 


^dr — Ad / dfc 


^DR — A d / dfc 


S(4S 2 -p 2 )’ 

- k D ~ 2 
t;(4t; 2 - p 2 ) 


E/(^) 


Note we need to perform the principal integration for to* 2 < p 2 /4. 
The pion decay constant at T = 0 reads the following simple form, 


f 2 _ N C M 0 

Jn DR — 


4 -D 




(2tt) d / 2 

We show a concrete example of trSj^ and / 2 DR for D ~ 2,3,4 in appendix [Bl 


(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 


(83) 


IV. MODEL PARAMETERS 

Having obtained the equation which determines the pion mass and pion decay constant, we are now ready to perform 
the parameter fitting. In the previous section we suppose that all the cutoff scales are equal in each regularization 
to reduce the parameters. Thus the model has three parameters: the cutoff scale A, the current quark mass m u and 
the coupling strength G. Whereas in the DR there appears one more parameter, so the total number becomes four: 
the current mass to u , dimension D : scale parameter Mq and the coupling G, as discussed in Ref. (36|. In this section, 
we shall set the model parameters by fitting the pion mass and decay constant. The actual values we use are shown 
below 


TOtt = 135 MeV, U = 94 MeV. (84) 

For the case with the DR, we perform fitting with one more quantity, the neutral pion decay constant to two photons, 
which will be discussed later. 


A. Parameters in various regularizations 

Here we align the model parameters in various regularizations in this subsection. 

Table n mi mu and m show how the parameters change according to the current quark mass m u , where we first 
set the value of m u then search the parameters A and G which lead m^ = 135MeV and /, r = 94MeV. One sees the 
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TABLE I. Parameters in the 3D cutoff 


m u 

(MeV) A 3D (MeV) G ■ 10” 6 (MeV 2 ) 

m*(MeV) (uu] 

) 1/3 (MeV) 

3.0 

942 2.00 

220 

-300 

4.0 

781 3.09 

255 

-272 

5.0 

665 4.71 

311 

-253 

5.5 

609 6.26 

375 

-245 



TABLE II. 

Parameters in 

the 4D cutoff 


m u (MeV) A 4D (MeV) G 

• 10" 6 (MeV“ 2 

) m* (MeV) (uu) 

1 ' /3 (MeV) 

3.0 

1397 

1.80 

198 

-300 

5.0 

1027 

3.64 

242 

-253 

7.1 

827 

6.66 

311 

-225 

8.0 

768 

8.88 

369 

-216 


TABLE III. Parameters in the Pauli-Villars regularization 


m u (MeV) Apv(MeV) G-10" 6 (MeV- 2 

! ) to* (MeV) ( 

mm) 1/3 (MeV) 

3.0 

1420 

1.77 

195 

-300 

5.0 

1071 

3.45 

229 

-253 

9.94 

780 

9.56 

311 

-199 

10.0 

778 

9.64 

312 

-198 

15.0 

729 

19.4 

417 

-173 


TABLE IV. Parameters in the proper-time regularization 

m u (MeV) A PT (MeV) G • Kr 6 (MeV- 2 ) 

to* (MeV) 

(mm) 1/3 (MeV) 

3.0 

1464 

1.61 

178 

-300 

5.0 

1097 

3.07 

204 

-253 

10.0 

755 

8.13 

265 

-198 

12.6 

680 

12.1 

311 

-183 

15.0 

645 

17.2 

372 

-173 


tendency that cutoff scale A becomes smaller with increasing m Ul while G becomes larger according to m u ■ We confirm 
that the values of the cutoff and four point coupling are 0(1) and O(10~ 6 ) in MeV scale in these regularizations. 

We also showed the values of the constituent quark mass to* and the chiral condensate (uu) 1 ' 3 which are the 
predicted quantities in the models. We note that the values of to* are about 200 — 300MeV which are comparable to 
one third of the proton mass. We also note that to* increases with respect to m u , while the absolute value of (uu) 
decreases when m u becomes large. Since the relation to* oc G(uu) holds, even (uu) becomes smaller m* can be larger 
due to the large value of G, which is actually the case in these regularizations. 

We plot the parameters of each regularization in Fig.l. The black circles denote the value which satisfy to* = 
311MeV for each regularization. The relation between the cutoff scale A and G for the 4D cutoff, Pauli-Villars 
regularization and proper-time regularization resembles each other. The relation between m u and A for these regu- 
latizations also resembles each other. In the case of same value for to*, the m u dependence of A is large. However, 
the values of A are close, 660 — 830MeV. In the case of same value for to*, the m u dependence on G is large and the 
values of G are separate in each regularization. The relation between m u and G is different in each regularization. 
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FIG. 1. Parameters in each regularization 


B. Parameter fitting in the dimensional regularization 


For the sake of the parameter fitting in the DR, we present the calculation of the pion to two photon decay rate, 
r 7 r o 77 , in this subsection. The decay rate can be evaluated through the following one-loop amplitude, ^2), 


f / ! V {h, k 2 ) = ^^ P \K k 2 ■ T i 


2 N c 


Try = g^qqe —m*M 0 


4 -D 


d u p 

i(2 tt) D (P — /C 2) 2 — TO* 2 p 2 — TO* 2 (P + fci) 2 — TO* 2 


(85) 

( 86 ) 


where e is the QED coupling constant and k\ and &2 are the external momentum of emitted photons so the square of 
the total momentum coincides with that of the original pion, namely, (fci + fe) 2 = to 2 . By using T 7 , the decay rate, 
r 7 r o 77 , is expressed as 


- 7T77 


-HhL\T 1 

" 64J 7 ' 


The detailed derivation is presented in the paper [27j. After some algebra, one obtains 


= -4 g. 


..N n - 


(4ttM 0 2 ) 


2\2-D/2 


T(2 — D/2) 


nqq 


37 t ( to * 2 ) 3_d / 2 


J dx - | [l - x(l - x)m^\ 2+D / 2 _ i| _ 


(87) 


( 88 ) 


with a e = e 2 /(47r), and to 2 = to 2 /to* 2 . 

With the observables, to v = 135MeV, = 94MeV and 


r w o 77 = 7.8 eV, (89) 

we perform the parameter fitting in the DR following [27] . Table. |V] shows the fitted parameters in the DR case. 
We note that both the constituent quark mass and chiral condensate grow up with increasing m u , which is the 
characteristic feature in this regularization |36] . 
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TABLE V. Parameters in the dimensional regularization 


m u (MeV) 

D 

GMq~ d ( MeV’ 

- 2 ) M 0 (MeV) 

m*(MeV) ( uu ) 

1/3 (MeV) 

3.0 

2.37 

-0.01134 

no 

-570 

-299 

5.0 

2.56 

-0.00588 

97 

-519 

-253 

8.0 

2.78 

-0.00241 

80 

-459 

-217 

20.9 

3.32 

-0.000229 

37 

-311 

-160 


V. MESON PROPERTIES 

We have presented the required equations in the model, then set the parameters for various regularizations. It 
is now ready for the actual numerical analysis on the model predictions. Here we shall show the thermal meson 
properties, which are the pion mass, the pion decay constant, and the sigma meson mass. 

At finite temperature, there are two ways of the application of each regularization; one is to apply the regularization 
only for the temperature independent contribution because the temperature dependent contributions are always finite 
due to the characteristic factor of the Fermi-Dirac distribution, i.e., f(E). The other is to apply the regularization 
both for the temperature independent and dependent parts, since the regularization essentially relates to the cutoff 
of the model so the introduction of the same cutoff clearly determines the model scale. On the other hand, the 
former method retains more symmetry of the model. Then the physical meaning of these prescriptions are that 
the former one respects the model symmetry, and the latter does the cutoff scale of the model. Since the model is 
not renormalizable, the predictions depend on the regularization ways and our purpose in this paper is to study the 
regularization dependence on the model. Therefore we shall study all the cases and compare the results among various 
regularizations. 


A. Results with regularizing T-independent contribution 

In this subsection, we show the results of meson properties based on the procedure of applying the regularization 
only to the temperature independent part. The required integrals are tr^ 1 and /(p 2 ), and we evaluate the following 
combination, 

tr S = trS^ eg + tiS T (90) 

where trSj( eg is trS for T = 0. The lower index indicates each regularization, namely, trS^, trS^j-,, trSp V , tr5p T and 
trS'{() R . For the temperature dependent part trS' T , here we use the form shown in Eq. EH) . Similarly for /(p 2 ), we 
use the equivalent expression, 


i(p 2 ) = iL s (p 2 ) + i t (p 2 ), 


(91) 


with I Reg for each regularization way. 

Figure 2 shows how the pion mass changes with respect to T and p for various parameter sets in the previous 
section. It should be noted that, for some parameter sets, no real solution exists at high temperature as seen in the 
case with the DR and m u = 3.0MeV [ 3 ll| . We observe the similar behavior in each regularization; the pion mass 
remains almost constant for low T and p, then raises up for higher T and p. This comes from the fact that the chiral 
symmetry is broken at low T and p and restores at high T and p. The pion has smaller mass when the symmetry is 
broken due to the Nambu-Goldstone theorem, while the mass becomes large after symmetry restoration. We see that 
the mass starts to increase around 170MeV which is comparable to the critical temperature for the chiral symmetry 
breaking. We see that the temperature and chemical potential where the pion mass glows up become larger with 
respect to m u for the 3D cutoff, 4D cutoff, PV and PT, while they become smaller for the DR. We also see that 
the discontinuity seen around the transition temperature is considerably larger in the DR case compare to the other 
regularizations. Then we expect that the tendency of the first order phase transition is strong for the DR comparing 
the other regularizations. We will discuss it in the next section. 

The results of the pion decay are shown in Fig. 3. One sees the similar tendency as well; the decay constant is almost 
constant for low T and p , and it decreases when T and p exceed certain values which are around T — 170MeV and 
p ~ 300 — 400MeV. It is interesting to note that the decay constant drops discontinuously at high p for some parameter 
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sets in the 3D, 4D, PV, PT regularizations, while the discontinuity is always the case in the DR. The existence of 
the gap is the signal of the first order phase transition, and the tendency becomes stronger with increasing m u . This 
is because the coupling strength is larger for higher m u as seen from the parameter tables, so quarks have stronger 
correlations when the parameter rn u is larger in the 3D, 4D, PV, PT cases. 

Having fixed the parameters with the pion mass and decay constant. We will calculate the sigma meson mass. It 
is one of the predictions of the model. Figure 4 shows the numerical results of the sigma mass. At T = 0 and fi = 0 
we find the band around 400 — 900MeV in 3D, 4D, PV, PT regularizations, 900 — 1400MeV in the DR. Then, in the 
DR case, the predicted values are larger than the experimental value nn a ~ 500MeV [37], in the leading order of the 
1/N C expansion. As is known we can obtain much smaller sigma meson width in this order. We should also check the 
next to leading order contributions for the sigma meson mass and width. The features of the curves are the similar 
to that of the pion; the mass decreases with increasing T and /i until some values, then it increases after exceeding 
the certain values. As seen in the pion mass case, the solution of the sigma mass on the real axis disappears for some 
parameter set at high temperature. 


B. Results with regularizing all contribution 


In this subsection, we study the meson properties by applying each regularization procedure to both temperature 
dependent and independent contributions. It is worth mentioning that the case with 4D cutoff method does not give 
credible results because enough number of frequency summations is not taken in this method. The cutoff scale of 
the 4D case is around lGeV, which means that at T = lOOMeV only four terms in the Matsubara mode summation, 
n = —2, ±1, 0 (u> n = (2 n + are taken into account. It is known that finite temperature field theory does not 

lead reliable predictions when the number of the frequency summation is small [38j]. Therefore, we will not show the 
results in the case of 4D cutoff scheme here, and consider the other four cases 3D, PV, PT and DR and call these cases 
3DRT, PVRT, PTRT and DRRT, respectively. It is also worth mentioning that the calculations technically become 
impossible at T = 0 in the PTRT case as can be read from Eqs. tEH) and m- Then, we will show the results with 
T = lOMeV as the representative values on n dependence at low T. 

Figure 5 displays the results of the pion mass. One sees that the qualitative feature does not change comparing 
to the previous case with regularizing only the temperature independent contributions. Quantitatively, we note that 
the changes become smoother at high T and /i. This can easily be understood because the regularization procedure 
suppresses the thermal contribution, so the finite temperature term reduces to give smoother curve with respect to T 
and fi. 

We aligned the results of the pion decay constant in Fig. 6. As seen in the pion mass case, the numerical results 
do not alter qualitatively, the curves become smooth. Note that, although the T dependence becomes considerably 
smother, the transition chemical potential does not change. This is due to the fact that finite temperature contributions 
become proportional to the step function 9(/j, — to*) for T = 0, then the transition chemical potential is not affected 
by the regularization procedure in this model treatment. 

The predictions on the sigma meson mass are shown in Fig. 7. The deviations from the results in the previous 
subsection in Fig. 4 are more or less similar to the deviations on the pion mass; the curves become flatter with 
respect to T while /i dependence does not indicate much difference for the 3D, PV and PT cases. However, there 
appears substantial difference between the cases of the DR and DRRT where the values of the sigma meson mass at 
fi = 500MeV are around 1500 — 2500MeV for the DR and around 600MeV for the DRRT case. This comes from the 
difference of the integral values between these two cases, which we will discuss in more detail in Sec. IVIIII 

We find that the solution on the real axis always disappears for high /r and low T(= lOMeV) in the PTRT with 
m u = 10, 15MeV. This is because for high /r some quantities become pure imaginary number in the calculation due 
to the complicated counter integral of /pt(.P 2 ) as seen in Eq. (1721) . Then one can not find a real solution in that case. 
This is the numerical reason why the meson properties behave badly for high /i in PTRT case. 


VI. PHASE DIAGRAM 


We shall draw the phase diagram in this section. We search the phase transition point by the condition that the 
maximum change of the chiral condensate with respect to T and [i. In more concrete, we numerically differentiate 
the condensate with respect to r e -^/T 2 + /r 2 , then the condition can be written 


d(uu ) 1 / 3 


= Max 


d(uu ) 1 / 3 


dr 


T—T c ,fj,=^, 


dr 


(92) 
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We should be careful on the case with the first order transition, because the condensate has the discontinuous point 
where we need to find the minimum of the thermodynamic potential then determine the chiral condensate. No matter 
whether the phase transition is the first order or cross over, we can also use the above criterion since at the first order 
point, d{uu) 1 ^ / dr = oo holds, therefore it is consistent in both cases. 

By searching the maximum number of the differentiate of the chiral condensate following the condition Eq. AMD, 
we draw a phase boundary for the chiral symmetry. Figure 8 shows the phase diagrams for the various parameter 
sets and regularizations. One sees that the critical temperatures, T c is found between 150 and 250MeV at g, = 0 
if we regularize only the temperature independent parts (left four panels and the bottom panel). On the other 
hand, a higher critical temperature is observed when the regularization is applied to the temperature dependent and 
independent parts. The regularizations 3DRT and PTRT give a critical temperature around T c ~ 200MeV at /r = 0, 
the PVRT induces a higher critical temperature near T c ~ 400MeV with m u = 15MeV, and the DRRT indicates it 
around T c = 300 — 400MeV. One also sees that the critical chemical potential, p, c , has no large dependence on the 
application of the regularization to finite temperature term, because the terms are dominated by the step function 
9(fx — m*) as mentioned in the previous section. Consequently, the area of phase boundary enlarges in the T direction 
when we apply the regularization to the temperature dependent term, which is numerically confirmed by the figure. 

We think the most interesting comparison from the figure is on the existence of the critical end point where the 
first order phase transition starts on the phase boundary. For 3D, 4D, PV, PT, no critical end point appears for the 
parameter sets with small rn u . While the critical end point always appears in DR. Then we can numerically conclude 
that the DR has stronger tendency of the first order phase transition comparing with the other regularizations. The 
existence of the critical end point in the other four regularizations can be understood by seeing the value of the 
parameter G. Briefly speaking, the critical end point appears when G is large. This is physically reasonable because 
G represents the strength of the correlation between quarks, then the larger G makes the condensation stronger. 


VII. PHASE DIAGRAM WITH FIXING m* 

We have seen the phase diagram and how the location of the critical end point depends on the parameters in the 
various regularizations. Considering the fact that the transition T and fi is essentially determined by the value of 
the constituent quark mass since its dependence appears in the thermal distribution, f(E), with E = y/k 2 + to* 2 , we 
think it may as well be interesting to compare the phase diagram with the parameter sets which lead the same value 
of m*(= 311MeV) at T = 0 and n = 0 instead of fixing the current quark mass, m u . 

First, we see the results based on the one with regularizing only the temperature independent parts, 3D, 4D, PV, 
PT and DR. The behaviors of m* for each regularization are shown in Fig. 9. All the results have no large difference, 
because the gap equations of the temperature independent part for each regularization has the similar behavior. The 
gap equations contain the following form, 

trS°(m*) 9 f(m*, A) + m* 2 ln(g(m*, A)), (93) 

in the regularizations, 3D, 4D, PV and PT. While in DR m* and A are replaced by Mq and D in some parts (see 
appendix iBl). In Fig. 9, we note that the results almost coincide in three regularizations, 3D, 4D and PT. The behavior 
in PV shows a smaller and DR gives steeper slope than the others. 

In Fig. 10 the phase diagram is illustrated in each regularization. The phase boundaries of 3D, 4D and PT show 
almost equivalent behavior. The area of the chiral symmetry broken phase for PV is larger than the others. This 
tendency comes from the behavior of m* in Fig. 9. Since the Apy is entered in the form of the dynamical mass, the 
chiral symmetry breaking contribution is enhanced than the other regularizations. The area for the chiral symmetry 
broken phase for the DR is smaller than the others. The critical end point for the DR locates higher temperature 
than the one for 3D, 4D and PT. These tendency also comes from the behavior of m* in Fig. 9. 

Next, we discuss the results with regularizing both the temperature independent and dependent parts, 3DRT, 
PVRT, PTRT and DRRT. The behavior of to* for each regularization is shown in Fig. 11. We find that the finite 
temperature effect becomes smaller or softer than 3D, PV, PT and DR, respectively. The behavior of to* in DRRT 
is the most closest behavior to the case of its regularizing only temperature independent part. 

The chemical potential has the similar contribution for 3DRT, PTRT, and PVRT. The behaviors of to* in DRRT 
and DR have large difference. This difference is caused by the reduction of the momentum integral dimension, D, 
from 4 to 3.32. 

The phase diagram in each regularization is shown in Fig. 12. To compare with Fig. 10, the region of the broken 
phase enlarges in PTRT, 3DRT and PVRT for a low chemical potential. However, the critical chemical potentials in 
PTRT, 3DRT and PVRT for a low temperature are almost equivalent to that in PT, 3D and PV, respectively. From 
the behavior of to* in Fig. 11, we observe a larger critical temperature and chemical potential for DRRT. We display 
the location of the critical end points (/tcpjTcp) for each regularization in Table IVTl 
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TABLE VI. Critical end point for each regularization 


Regularization fj, cp (MeV) Tcp (MeV) 


3D 

330 

25.0 

4D 

333 

46.8 

PT 

330 

26.0 

DR 

289 

74.7 

3DRT 

330 

25.0 

PTRT 

332 

23.2 


VIII. DISCUSSIONS 

We have introduced the regularization methods then studied the meson properties and phase diagram in previous 
sections. In this section, we are going to present more detailed discussions on the obtained results. 

In comparing the left panels v.s. right panels in Fig. 8, one can observe the contribution to apply the regularization 
procedure to the thermal correction. The phase diagram does not show considerable difference for 3D and PT cases, 
while in PV the area of the phase boundary becomes larger in the right panel when mn u is large, and the areas in 
DRRT case is larger than DR case for all the parameter sets. These difference can be understood through the following 
discussion. 

The loop integral, I, essentially has the following subtracted forms for 3DRT, PVRT, and PTRT, 


I 3 DRT = 

pOO poo 

/ d 3 fc F(k) — / d 3 k F(k), 

(94) 


Jo 

J A 


2pvRT = 

/ d 3 k [F{k,m*) - F(k, A)], 

Jo 

(95) 

2fpTRT = 

poo 

/ d 3 k 
Jo 

F[k) - F{k)( 1 - e -^ 2 +m- 2 )/A 2 j 

(96) 


Where the typical form of F(k,m*) is given by F(k,m*) = Cf(E)/E with some constant value, C. Thus the 
subtracted terms basically relate to the suppression on the high energy contributions which are expected to be small. 
We numerically confirmed that the subtracted parts are small for almost all the cases, then the phase diagrams does 
not change drastically. However, in the PV case with large m u , the difference between F(k, m*) and F(k , A) is small 
since the constituent quark mass becomes comparable to the cutoff scale, e.g., m* = 417MeV and Apy = 729MeV for 
m u = 15MeV in PV. Consequently, the thermal contribution strongly suppressed in the PV case with large m u . 

We saw that the area of the phase boundary does not alter so much in 3D, PV and PT regularizations. The essential 
reason is that the infinities appearing from loop integrals are subtracted at high energy. However in DR, the situation 
is different since the integral is replaced as 


/ d 4 k F(k) -A / d D k F(k), (97) 

Jo Jo 

so this modifies the integral kernel rather than the subtraction of high energy modes. This is the reason why DR shows 
considerable difference, if we apply the regularization to both temperature independent and dependent contributions. 

We also saw the location (or existence) of the critical end point is non-trivial. In Fig. 10, the diagram has the 
critical end point for 3D, 4D, PT and DR. No critical end point appears in the PV. Thus we find that the PV has 
weaker tendency of the first order phase transition than that of the other regularization methods. Particularly, the 
temperature of the critical end point in the DR case is higher then others, which may enable us to conclude that the 
DR has the stronger tendency of the first order phase transition. 


IX. CONCLUDING REMARKS 

We have studied the regularization dependence on the phase diagram of quark matter on T — /i plane by using the 
NJL model. We have first presented the regularization procedure at finite temperature and chemical potential, then 
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fitted parameters within various regularization methods. Thereafter, we have studied the meson properties and the 
phase structure. 

We find that the model produces the reliable predictions on the meson properties whose behavior for finite temper¬ 
ature and chemical potential does not alter drastically, which indicates that all the regularizations employed in this 
paper nicely capture physics on the meson properties. We can conclude that the regularizations are safely adopted. 
In this context the regularization parameter independent approach is also interesting [39,03] • 

It is expected that observation of the critical end point can distinguish a suitable regularization for an effective 
model of QCD by comparing with the resulting phase diagrams. The important difference is the existence of the 
critical end point. The model predicts that the critical end point appears at intermediate chemical potential around 
H ~ 300 — 400MeV. This density coincides the one in which different quark state such as color superconductivity 
may occur, there the order of the phase transition might affect crucially on such the dense states. Moreover the color 
superconductivity may be realized in the dense stellar objects, like quark stars and neutron stars jUJ. Therefore the 
study of the order of the phase transition has important meaning as well in cosmological observations. So we believe 
that the further and more extensive investigations are necessary on this subject. 
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Appendix A: Analytic expressions for I°(p 2 


Since the I°(p 2 ) integral can be evaluated analytically, we will present the explicit expression for various regular¬ 
izations. 

One needs special care in performing I(p 2 ) integral since it contains divergent contribution as seen in the 3D cutoff 
scheme. I® D becomes for m* 2 > p 2 / 4, 


1° (A 2 ) - 

) ~ 47r 2 


ln^D. 


*2 


m”“ , / 1 

+ 4a arctan | — 
2 a 


— 4b arctan ( — ] — 


2A? 


4D 


ap z 


arctan ( — 
2 a 


and for m* 2 < p 2 / 4, 


jo („2\ _ 

) - 47r 2 


ln(A| D + m* 2 ) + 4 a arctan ( — 

\ 2 a 

+ 2 - 2 ^fb ±a V n fb ±a 


where 


a = 


' A 4E> 


m* 


b = 


I A 4D 


Concerning on Ip V (p 2 ), it is convenient that we divide the integral as 

tO rO(m) .0(A) 

4py — J PV J PV > 

where Ip ( y FV ' > is subtracted part in the original integral and it becomes 


.0(A) 

'PV 


N, 


_ Jv c 

27T 2 


± L 


- ± c In - ± c 


with 


c = 


(Ap V — to* 2 + k 2 ) 2 A| 


PV 


4p 4 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 


(A6) 
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As seen above, we need to separately evaluate the integral Jpy" ' 1 
for to* 2 3 4 5 6 7 8 9 10 11 12 > p 2 /4, 


depending on the values of to* 2 and p 2 . It becomes 


jO(m) _ 

j pv ~ 


Nc_ 

47T 2 


In (to* 2 ) 


4 d arctan 



and for to* 2 < p 2 / 4, 


where 



(A7) 


(A8) 


(A9) 


Appendix B: trSn R and f n for D ~ 2,3,4 


We arrange the concrete expressions for trSpp and 
For D ~ 2, {D = 2 + 2e) we have 


- ^2 m * ■ M " 


1 171 *^ 

-+ lE -]n(2n) + \n ^ 


IV, 




1 + e < 7s — ln(27r) + In 


TO 


*2 


For D ~ 3, (D = 3 + 2e) we have 


trSoR — ^- C to* • \/ m* 2 M 0 


fn dr — 

For D ~ 4, (D = 4 — 2e) we have 


y/2w 

N c 

2V2* 


2 TO* 2 

1 + e ^ 7 B + In — — 2 + In —j 


V m* 2 Mn 


2 777,* 2 

l + e{ 7 B + ln- + l„— 5 - 




Af c 

47T 2 ’ 

IW 


"* 2 - - 7 s + ln(27r) + 1 + In 

e m z 


^DR- 47r 2 TO 


M 2 


-7e + ln(27r) + In —» 


m 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 
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FIG. 2. Pion mass. Left: fj, = 0. Right: T = 0. 
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FIG. 3. Pion decay constant. Left: fi = 0. Right: T = 0. 
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FIG. 4. Sigma mass. Left: /j, = 0. Right: T = 0. 
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FIG. 5. Pion mass. Left: /x = 0. Right: T = 0 for 3DRT, PVRT, DRRT and T = lOMeV for PTRT. 
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FIG. 6. Pion decay constant. Left: n — 0. Right: T = 0 for 3DRT, PVRT, DRRT and T = lOMeV for PTRT. 
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FIG. 7. Sigma mass. Left: n = 0. Right: T = 0 for 3DRT, PVRT, DRRT and T = lOMeV for PTRT. 
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FIG. 8. Phase diagram. 
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FIG. 9. Comparison of the constituent quark mass with fixed parameters under m* = 331MeV. Left: 


fj, = 0. Right: T = 0. 



FIG. 10. Comparison of the phase diagrams with fixed parameters under to* = 331MeV. 




FIG. 11. Comparison of the constituent quark mass with fixed parameters under m* = 331MeV. Left: fi = 0. Right: T = 0 
for 3DRT, PVRT, DRRT and T = lOMeV for PTRT. 
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FIG. 12. Comparison of the phase diagrams with fixed parameters under m* = 331MeV. 






